Full separability criterion for tripartite quantum systems 
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In this paper, an intuitive approach is employed to generalize the full separability criterion of 
tripartite quantum states of qubits to the higher-dimensional systems (Phys. Rev. A 72, 022333 
(2005)). A distinct characteristic of the present generalization is that less restrictive conditions are 
needed to characterize the properties of full separability. Furthermore, the formulation for pure 
states can be conveniently extended to the case of mixed states by utilizing the kronecker product 
approximate technique. As applications, we give the analytic approximation of the criterion for 
weakly mixed tripartite quantum states and investigate the full separability of some weakly mixed 
states. 
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I. INTRODUCTION 



Entanglement, as an essential ingredient of quantum 
information theory, has been an important physical re- 
source for a lot of quantum protocols, such as quantum 
computation [1], quantum cryptography [2], quantum 
teleportation [3], quantum dense coding [4] and so on. 
Recently, many efforts have been made to characterize 
the quantatively properties of entanglement [5-8], how- 
ever, the good understanding is only restricted to low- 
dimensional systems. The quantification of entanglement 
for higher dimensional systems and multipartite quantum 
systems remains an open question. 

Since Coffman et al. [9] introduced the so called resid- 
ual entanglement on the basis of concurrence [5], the 
investigation of multipartite entanglement has attracted 
much attention. For example, Diir et al. have considered 
the classification of entanglement for tripartite systems 
of qubits [10]; Miyake [11] has given the classification 
for multipartite systems based on the hyperdeterminant. 
On the basis of the different classes of multipartite entan- 
glement, the corresponding entanglement monotones can 
be given [11,12]. Some quantities have also presented to 
characterize the properties of entanglement by collecting 
the contributions of the entanglements of different classes 
[13,14]. One can note that the quantities introduced in 
Refs. [13,14] can also characterize the full separability 
of a pure multipartite state. However it is easily found 
that construction of these quantities requires more re- 
strictive conditions. Even though some conditions may 
be repeated, it is usually not easy to exclude the repeated 
ones, especially for high-dimensional systems. Hence, it 
will reduce the efficiency of calculation to some extent. 

Considering the full separability criterion introduced 
in Ref. [15], which can effectively reduce the restrictive 
conditions to some extent, in this paper, we will gen- 



eralize the criterion to high-dimensional systems by an 
intuitive approach. The generalized full separability cri- 
terion for pure states can be conveniently extended to the 
case of mixed states by utilizing the kronecker product 
approximate technique which can usually further reduce 
restrictive conditions. As applications, we give the an- 
alytic approximation of the criterion for weakly mixed 
tripartite quantum states and study the full separabil- 
ity of some weakly mixed states. The paper is organized 
as follows. Firstly, we give the intuitive generalization 
of the separability criterion for pure states; secondly, we 
extend it to mixed states and discuss the full separability 
of some quasi pure states; the conclusions are drawn in 
the end. 



II. FULL SEPARABILITY CRITERION FOR 
TRIPARTITE PURE STATES 

At first, let us recall the full separability criterion for 
tripartite pure states of qubits given in Ref. [15]. A 
tripartite pure state {tp) abc denoted by a vector in 2 x 
2x2 dimensional Hilbert space, 

IV') = (flOOOi flOOlj iQlOj iQllj Q-lOOj Q-lOlj aiio,aiii)"^j 

with the superscript T denoting transpose, is fully sepa- 
rable, if and only if 
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C(|V')) = |C(|V^))| = ^X^|C"|^ = 0, (1) 

here the vector C(-0) = © C" with C" = (i/'*|s" [i/-), 
where the star denotes complex conjugation, and 

— —ay®(Ty®Ii, = -ay®ay®l2, — ~ay®Ii®ay, 

(2) 

S'^ = -(Jy®l2®(Ty,S^ = -h®ay^ay,S^ = -h^dy^CTy, 

(3) 
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As mentioned in Ref. [15], a tripartite pure state of 
qubits ean be considered as a tensor cubes. Directly, 
a tripartite higher- dimensional pure state can naturally 
considered as a tensor grid which includes tensor cubes. 
E.g. let|0^Bc)=EL=oEfe =0 ^ijk \ijk)jiBc^ tensor 
grid of \(j)ABc) can be pictured as two adjoining cubes, 
which includes three tensor cubes. In this sense, one can 
draw a conclusion that tensor cube can be regarded as the 
unit of tensor grid. Since every tensor cube in a tensor 
grid can be considered as an non-normalized tripartite 
pure state of qubits, one can get that every unit corre- 
sponds to a C defined in eq. (1). Therefore, the tensor 
cube can also be considered as a unit which describes the 
full separability of a tripartite higher-dimensional pure 
state. In other words, the full separability of the given 
tripartite higher-dimensional pure state can be described 
by the full separability of the non-normalized tripartite 
pure state of qubits. 

Theorem l:-For any a tripartite pure state |x) which 
includes M non-normalized tripartite pure states of 
qubits (tensor cubes mentioned above), let the non- 
normalized pure state of qubits corresponding to the ith 
cube be denoted by \(pi), one can obtain the correspond- 
ing C (|<^i)). Define 
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M 



(5) 



for the state |x), then |x) is fully separable, if and only 
ifF(|x)) = 0. 

Proof. It is obvious that C(|x)) = means that 
C' = holds for all ipi, vice versa. Since the tensor 

cube corresponds to the unit of describing full separabil- 
ity, C(|x)) = shows that there does not exist any en- 
tanglement in Ix). That is to say, the tripartite quantum 
state Ix) is fully separable. In other words, since every 
non-normalized (pi is fully separable, one can obtain that 
every group of parallel lines of the tensor grid is linear 
dependent. I.e. the state that the grid denotes is fully 
separable [15]. On the contrary, if |x) is fully separable, 
C{\^i)) = 0, i.e. C(|x)) = 0. 

Considering the matrix notation of 

ni — 1 712 — 1 n^ — l 

\x) = J2 12 1*-^'^^ ' 

i=0 j=0 fe=0 

C(|x)) can be expressed as the function of |x), i-e. 



cm = 



Ni N2 Ns 



EEE^'((«a®«/3®«7)lx)), (6) 



where Np = with p = 1,2,3; Sq, q = a,/3,7, 

denotes 2 x Up matrix with p corresponding to q. If the 
generator of the group SO{np) is denoted by (S'p, Sq can 
be derived from |S'p| by deleting the row where all the 
elements are zero, where | | denotes the absolute value of 
the matrix elements. 

According to eq. (1), eq. (6) can be expanded by 



C{\X)) = 



Ni N2 N3 9 

\ EE EE i^'((««®^/5®«7)ix))i' 

\ a=l /3=l7=l 5=1 

Ni N2 N3 9 

EEEE((x15j^,S^^a/57lx) 
a=l /3=1 7=1 5=1 

1/2 



X {x\S^f,^s'S^0^\x*))] 



(7) 



where Sa/s-y = Sa ® s/3 ®) s^, s^ are defined by eqs. (2-4), 
and the superscript T denotes transposition operation. 



III. FULL SEPARABILITY CRITERION FOR 
MIXED STATES 

On the basis of C(|x)) for pure states, the corre- 
sponding quantity C{p) for mixed states p defined in 
Cdxd{d = ni X n2 X nz) is then given as the convex of 



C(p)=inf^ftC(|*,)) 



(8) 



of all possible decompositions into pure states |\l/i) with 



P = ^Pi\^i) >o. 



(9) 



C{p) vanishes if and only if p is fully separable. Substitute 
eq. (7) into eq. (8), one can get 



C{P) 



inf^p. 



Ni N2 N3 9 

EEEEI^'(^«/^7i*.))| 

a=l /3=1 7=1 S=l 



1 1/2 



a=l /3=l7=l 



(10) 

It is obvious that if the infimum of eq. (10) can be pro- 
vided, one can obtain a sufficient and necessary condition 
of separability for mixed states. However, it seems to be 
impossible for higher-dimensional systems. One can only 
give a lower bound as a necessary condition. Therefore, 
a lower bound with strong sufficiency or convenient for 
calculations is expected. 

According to the matrix notation [7] of equation (9), 
one can obtain p = ^W'^^, where W is a diagonal matrix 
with Wii = Pi, the columns of the matrix ^' correspond 
to the vectors |\l/i). Due to the eigenvalue decomposition: 
p = $M$'^, where M is a diagonal matrix whose diagonal 
elements are the eigenvalues of p, and <I> is a unitary 
matrix whose columns are the eigenvectors of p, one can 
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obtain l-VF^/^ = ^Afi/^JJ, where [/ G C""'^^ is a Right- 
unitary matrix, with N and r' being the column number 
of ^ and the rank of p. Therefore, based on the matrix 
notation, eq. (10) can be rewritten as 

C{p) ^ inf./ V |C/^MV2$T5T 55^^^$^i/2f/' 



where the Minkowski inequality 



\ 



Pm^n 



is used. According to Ref. [7], one can di- 
rectly obtain a lower bound of C{p) as max Ai(z) — 

z 

J2i>i^ii^)' where Xj{z) are the singular values of 
Ei=i ZjM'^/'^^'^S'^i^^s^Sa^j^M'^/'^ in decreasing order 
with z = [zi, Z2, • • •, z/] a group of optimal complex pa- 
rameters. It can be easily found that the number of op- 
timal parameters (/ = 9iVi ■ N2 ■ N3) is too large to be 
conveniently used to calculations for higher-dimensional 
systems yet. However, it will be found that by kronecker 
product approximation technique, not only might the 
number of optimal parameters be further reduced, but 
also one can calculate the lower bound in different ap- 
proximation degrees. In particular, we can provide an 
analytic approximation for weakly mixed states. 

In fact, if replacing "x" of eq. (11) by "O", eq. (11) 
can be rewritten as 



Cip) > inf{ Y: [{U^M'/H^SI,^S'S^,,^M'/'U) 

0!.,f3.'y.5 



m{{[{U^^U^)A{U^U*)y;r}y\ (12) 



where 



^ = E E (P^/^)^ Sj,,E^S„,, (pV2) , (13) 



a,/3,7 5—1 



defined in CdxdOQxd, and pi/2 = (<^M^/^)(^{<^M^/^)* , 

Sa/97 = Safj^f (S Safjj- The othcr indices in 
above equation arc all defined the same as previous sec- 
tions. Even though the value of eq. (11) is not changed, 
the implied meaning is quite different, which means that 
we have copied the given quantum state in a conjugate 
Hilbert space and we consider the separability of the state 
in a doubled Hilbert space. The distinct advantage is 
that eq. (12) allows us to employ the kronecker product 
approximation technique [16,17]. 



Next we will employ the kronecker product approxima- 
tion technique on A to derive a lower bound of eq. (12). 
Based on the technique, A should be converted [19] into 
iby 

A = Vi2{AVi2f', 

where the superscript T2 denotes partial transposition 
on the second subspace [18], V12 is swap operator [19] 
defined as 

ikj' k' 

j,k' = 1,- • ■,d,j',k = 1,- • -,(1. 
A has the singular value decomposition: 



(14) 



i=l 



where Ui, Vi are the ith columns of the unitary matrices U 
and V ^ respectively; S is a diagonal matrix with elements 
(7j decreasing for i = 1, • • •, r; r is the rank of A. Thus, 
based on Ref. [16,17] A can always be written by 



i i 

where Vec{Ai) — y/ffiUi and Vec{Bi) = ^/cTiv*. For any 
a.px q matrix M = [niij] with entries [20], Vec{M) 
is defined by 

Vec{M) = [mil, •••,mpi,TOi2, •• •, mp2, • • •, Wig, mp,]^. 

(15) 

One can find from eq. (13) that if the two subspace that 
A is defined in is exchanged, A will converted into A*, 
hence one has Bi = , i.e. 



(16) 



Substitute eq. (16) into eq. (12), eq. (12) can be given 
by 



Af r 



c(p)^i?;EEl(t^^A-^)i 



(17) 



The infimum can be employed to test the full separability 
of p. 

In terms of the Cauchy-Schwarz inequal- 

1/2 / ^ 1/2 



^ Y^^iVi and 



Eikil > Ei^Jij, C{p) given by eq. (16) can ar- 
rive at 



N 



C(p)>infE 



i=i 



(18) 
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where Zj = Xj cx:p{i4>j), with Xj > 0, '^j^'j = 
1. Therefore the infimum of eq. (18) can be given by 
max Ai(z) — X^j>iAj(z), where Xj{z) are the singu- 
lar values of ^X]j=i ^j^j^ decreasing order [7], with 

z — [zi,Z2,- ■ -tZr]- Note that r < is usually much 
smaller than cP in practical calculations. In particular, 
one can consider different numbers of aj in decreasing or- 
der and correspondingly introduce optimal parameters, 
which will might provide approximate lower bounds in 
different degrees. In this sense, the number of optimal 
parameters can be dramatically reduced. In fact, it is 
very possible that Aj corresponding to the maximal aj 
can give the main contribution [8] to the infimum of eq. 
(18). That is to say the lower bound of C{p) can be given 
by Ai — Aj with A^ the singular values of Aj, which 
is an analytic approximation. 

For weakly mixed states i.e. quasi pure states, an an- 
alytic approximation of C(p) can also be introduced [21]. 
According to eq. (13) and kronecker approximation tech- 
nique, A can also be given in the following way 

9 

a,/3,7 (5=1 

X ((*?! sl^^s's^p^ X slp^s's^p^ I*;,,)) , 

(19) 

where ^'q, and Ua denote the ath eigenvector and eigen- 
value, and all the other quantities are defined similar to 
those in eq. (7). According to the symmetry of A given 
by eq. (16) and the kronecker product approximation 
technique in above section, A can be formally written as 

Aim \ ^ rpot (rnoi \* 

■^I'm' — 2^ ^ Im l-* I'm' ) ■ 
a 

The density matrix of quasi pure states has one single 
eigenvalue /zi that is much larger than all the others, 
which induces a natural order in terms of the small eigen- 
values > 1. Due to the same reasons to those in Ref. 
[21], here we consider the second order elements of type 
j4'™. Therefore, one can have the approximation 

A"^' — TlmTiim' ^ith Tim = ~~Ffff- 

V ^11 

In this sense, eq. (18) can be simplified significantly: 

C(p)c^Ca(p) = inf^|C/^TC/|... 

i 

Ca{p) can be given by 

Ca{p) = max{Ai - ^ Ai,0}, 

i>l 

where A, is the singular value of r in decreasing order. 



Consider two (2x2x3) —dimensional quasi pure 
states constructed respectively by 

pi{x) = X \GHZ') {GHZ'\ + (1 - x)li2 

and 

P2{x)=x\W') {W'\ + {1-X)U2, 

where 

\GHZ') = i(|000) + jlOl) + [Oil) + |112)), 

and 

|I^') = ^(|000) + |011) + |112)). 

Note that \GHZ') and \ W') given in Ref. [11] correspond 
to GHZ class and W class with high local rank, respec- 
tively. The two states can be considered as quasi pure 
states for x > 0.3. By the calculation, one can find that 
Co(pi) and Ca{p2) are both nonzero. What is more, for 
the quasi pure states generated by the mixture of maxi- 
mally mixed state (identity matrix) and tripartite GHZ 
state in 3 X 3 X 3 dimension, the corresponding Ca{p)s can 
all be shown to be nonzero for x > 0.3. We also study 
some (2x2x3) —dimensional quasi pure states p by the 
mixture of maximally mixed state and random semisepa- 
rable pure states generated by Matlab, numerical results 
show that Ca{p) are nonzero if p are strict quasi pure 
states. All above show the sufficiency of our criterion 
for testing the entanglement of high- dimensional mixed 
systems. 

IV. CONCLUSION AND DISCUSSION 

In summary, we have utilized an intuitive approach 
to generalize the criterion to high-dimensional tripartite 
systems. The generalized criterion for pure states can 
be conveniently extended to the case of mixed states by 
utilizing the kronecker product approximate technique. 
The lower bound for mixed states can provide necessary 
conditions to test the full separability. Compared with 
the previous criteria, the criterion introduced here can ef- 
fectively reduce the restrictive conditions. However, the 
criterion is not an entanglement monotone. Numerical 
results show that our criterion for high-dimensional sys- 
tems is even sufficient condition of full separability for 
strict quasi pure states. 
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